Narrow-escape-time problem: the imperfect trapping case 
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We present a master equation approach to the narrow escape time (NET) problem, i.e. the 
time needed for a particle contained in a confining domain with a single narrow opening, to exit 
the domain for the first time. We introduce a finite transition probability, u, at the narrow escape 
window allowing the study of the imperfect trapping case. Ranging from to do, v allowed the study 
of both extremes of the trapping process: that of a highly deficient capture, and situations where 
escape is certain ("perfect trapping" case). We have obtained analytic results for the basic quantity 
studied in the NET problem, the mean escape time (MET), and we have studied its dependence 
in terms of the transition (desorption) probability over (from) the surface boundary, the confining 
domain dimensions, and the finite transition probability at the escape window. Particularly we show 
that the existence of a global minimum in the NET depends on the 'imperfection' of the trapping 
process. In addition to our analytical approach, we have implemented Monte Carlo simulations, 
finding excellent agreement between the theoretical results and simulations. 

PACS numbers: 05.40.Fb 



I. INTRODUCTION 

The time needed for a particle contained in a confin- 
ing domain with a single small opening to exit the do- 
main for the first time, usually referred as narrow es- 
cape time problem (NET), finds a prominent place in 
many domains and fields. For instance in cellular biol- 
ogy, it is related to the random time needed by a particle 
released inside a cell to activate a given mechanism on 
the cell membrane ([H-Q)- Generally speaking the NET 
problem is part of the so called Intermittent processes, 
which are used to explain scenarios ranging from animal 
search patterns (0|), through the solutions or melts of 
synthetic macromolecules ([|, @), to the manufacture of 
self-assembled mono- and multi-layers (@, H|). Since the 
work of Berg and Purcell (|9|), research in the NET prob- 
lem area has experienced a stea dy g rowth over time and 
motivated a great deal of work ( [T0l - [2ll |). 

In [r|, we have introduced an analytical Markovian 
model that showed the impact of geometrical parameters 
and the interplay between surface and boundary paths in 
the studied confining domain, of a discrete and rectan- 
gular shaped nature, for the perfect trapping case. With 
"perfect trapping" we refer to the particle's impossibility 
of return to the system, i.e. once the particle reaches 
the narrow opening the escape becomes certain. In that 
work we presented a phase diagram which showed that 
some combinations of the geometrical parameter and the 
transport mechanism were required for the existence of 
an optimal transport (a global minimum in the NET). 

In this work we consider the same confining domain 
and the transport properties that we dealt with in Ref. 
[Tf^ |. However we eliminate the assumption of perfect 
escape introducing a finite transition probability at the 
narrow escape window. It is well known that systems de- 
scription through the "imperfect trapping case" (once in 
the trap site capture is not certain) are suitable whenever 
the surface contains 'deep traps', capture and re-emission 



from a surface that contains sites with several internal 
states such as the 'ladder trapping model', proteins with 
active sites deep inside its matrix, etc. (|22h25|). Under 
this new assumption, we have discovered some very inter- 
esting results. Particularly we show that the existence of 
a global minimum in the NET depends on the existence 
of an imperfection in the trapping process. 

The aim of this work is to study the influence of the 
'imperfection' in the passage through the escape window 
on the effective diffusion process, and its effect on the 
NET problem. For that purpose we calculate the time re- 
quired by the particle to escape from the system (MET) . 
In order to perform our research we exploit Dyson's the- 
ory ([HI]), and the notions of Absorption Probability Den- 
sity (APD) (|1|). 

The outline of this paper is as follows. In the follow- 
ing section we introduce some general results regarding 
imperfect trapping process as well as our model, and pro- 
vide the basic definitions and concepts. We also describe 
the proposed analytical approach and present the main 
results. Section III depicts several assorted illustrations 
for the MET to the narrow opening for different configu- 
rations of the system through a comparison between our 
analytical framework and Monte Carlo simulations. In 
Section IV we discuss our conclusions and perspectives. 
Finally, in the Appendix we present the calculations in- 
dicated in Section II. 



II. ANALYTICAL APPROACH 

A. Some general results regarding imperfect 
trapping 

Let us consider the problem of a walker making a ran- 
dom walk in some finite domain with a trap or sink 
present in the system. We will follow the walker evo- 
lution through the system considering the 'unrestricted' 
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conditional probability P(s, t\so, t = 0), that is, the prob- 
ability that a walker is at s at time t given it was at Sq 
at t = 0. By 'unrestricted' we identify a situation with 
no traps/sinks present in the system. 



V 



FIG. 1. Finite domain with a trap or sink present. The 
entrance to the trap site (empty circle) is regulated at the 
escape 'window' by the transition rate v. 



Absorption Probability Density and Mean Absorption Time 

As we are interested in the trapping process, let us de- 
fine A(s,t\so,0) as the Absorption (trapping) Probability 
Density (APD) through the site s at time t, given that 
the walker was at sq at time t = 0, i.e., A(s, t\so, 0) dt 
gives the trapping probability of the walker, through s, 
between t and t + dt given that it started at t — from 
sq. It is worth mentioning that the First Passage time 
approach is not fully applicable since an 'excursion' to 
the trapping site does not necessarily ends the process. 
However we show in the following lines that an inter- 
esting relation could be established between A(s, t|s*o,0) 
and F(s, t\s , 0), the First Passage time density (FPTD) 
through the site s at time t, given that the walker was at 
s*o at time t = 0. 

From now on we will denote the Laplace transform of 
a function f(t) by its argument, 



£{/(*)} = /(«) 



The connection between the APD and the 'unre- 
stricted' conditional probability P(s, t\sp,t = 0) could 
be traced to results in [23| or [27]] • This approach in the 
Laplace domain gives, 



A(0,u|s* 0) t = 0) = 



vP(0,u\s o ,t = 0) 

i+vP(d,u\d,t=o) 



(i) 



Let us make a brief digression regarding the relation 
between the APD and the FPTD. For this we rewrite Eq. 
([1]) in the form 



A(0,u\s o ,t = 0> 



P(6,u\s ,t = 0) 

" p(6,u\6,t = o) v + 



(2) 



P(0,«|0,t=0) 



As usual the connection between FPTD and the 'unre- 
stricted' conditional probability P(s,t\so,t = 0) is es- 
tablished (in the Laplace domain) through the 'Renewal 



approach' ([28]), 



F(6,u\s ,t = Q) = 



P(6,u\s ,t = 0) 



(3) 



P(0,u\0,t = 0) 

We recognize in equation ([2]) , the FPTD ((3]) and using the 
relation F(6,u\6,t = 0) = 1 - *(0,u)P(0,w|0,t = 0) _1 , 
where ^(s, r) is the probability that the walker remains 
at s until time r since it arrived at s at time (in the 
unrestricted case) [2^ |. we rewrite Eq. © as 

A(0,u\s o ,t = 0) = F(6,u\s ,t = 0) 
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*(0,m)- 1 +i^ 



F(d,u\6,t = 0) 



(4) 



The term ^(^(0, u) 1 + v) 1 in ((4} gives the fraction of 
walkers that are trapped while ^(0, u)(^(0, u)^ 1 + u)^ 1 
gives the ones that are not absorbed. Further considera- 
tions could be made regarding (U]), we write it as, 



A(0,u\s o ,t = 0) = ^2A-j(6,u\s 0l t = 0) 

3=1 



where 



A J (0,u\s o ,0) = F{0,u\s o ,0)\ 



f{o,u\o, o)^(o, u y 



(5) 



3-1 



(6) 



Notice that Aj(-) accounts for that walkers that are ab- 
sorbed in the j-visit (j = 2,3,...) and not before, i.e., 
the walkers arrive for the first time at site from sq but 
these are not absorbed until they return to site for the 
(j — l)-time. 

The probability of being absorbed in the j— visit at site 
is, 



Aj(d,t\s ,0)dt= Aj(6,u = 0\s ,0) 



(7) 



^(0,0)-! + 



3-1 



where ^(0, u = 0) is the mean residence time at site in 
the unrestricted system and we have used that F(s, u — 
0|s*o,t = 0) = 1 for a finite (unrestricted) system. As 
equation (O shows, when v — > independently of the j 
value, we have no absorption, while in the limit v — > oo 
the absorption is certain in the first 'visit' to the site (i.e. 
3 = !)• 

The mean absorption time, i.e., the mean time until 
the walker is absorbed is evaluated in terms of A(-) as, 

T= / t J2A(d,t\s ,0)g(s Q )dt 
Jo - 



^A(0,m|s ,0) 5 (s ) 



(8) 



u=0 
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where .g(sn) denotes the probability density of initially 
finding the walker at a position sq ([28]). 



B. The Model 

For our model we consider the problem of a walker 
making a random walk in a finite rectangular N x (M+l) 
lattice (see figure [2]). The surface is bounded in the y di- 
rection where the walkers can move from y — to y = M, 
and periodic boundary conditions are assumed in the x 
direction so x and x + N denote the same place in space. 
As we mentioned before, we follow the walker's evolution 
through the system considering the conditional probabil- 
ity P(n,m,t\no,mo,t — 0) = P(n,m,t), where (n,m) 
are discrete coordinates in the {x,y) space. P(n,m,t) 
satisfies the following master equation: 

P(n, 0, t) = jP(n, 1, t) - 5P{n, 0, t) 

+/3(P(n + 1, 0, t) + P(n - 1, 0, 0, t) 
-2P(ra, 0, f)) ; m = 

P{n, 1, t) = SP(n, 0, t) - 4 7 P(n, 1, t) 

+ 7 (P(n +l,l,i) + P(rc- 1,1,*) 
+P(ra,2,t)); m = l 

P(n, m, t) = j(P(n - 1, m, t) + P(n + 1, m, t) 
+P(n, m + 1, t) + P(n, m - 1, i)) 
-47P(n, m, t) ; 2 < m < M - 1 

P(n, M, t) = j(P(n - 1, M, t) + P(n + 1, M, t) 
+P(n,M -l,t)) 
-3jP(n, M, t); m = M 

(9) 

where 7 is the surface transition probability per unit time 
in the x and y direction, /3 is the transition probability 
over the line m = in the x direction, and <5 is the 
desorption probability per unit time from the boundary 
line m = 0. 

We introduce the imperfect escape case by allowing 
a finite transition probability {y) at the narrow escape 
window. Varying from to 00, v allowed the study of 
both a deficient trapping, and situations where escape is 
certain (i.e. perfect trapping case). 

In the following we will say that the walker 'escapes' 
when it gets trapped or adsorbed, without the possibility 
of returning to the system. This terminology matches the 
one used in the NET area. Hence, adsorption — > escape, 
and so on. 



Escape Probability Density (EPD) 

We now make a brief comment regarding the Escape 
Probability Density. Taking into account the parameters 
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FIG. 2. Schematic transitions of the walker to/from the base 
line and to/from a generic surface site. Notice that the en- 
trance to the trap/escape site (empty circle) is regulated by 
the transition rate v also notice that it could be reached both 
from the surface (with transition rate 7) and from the baseline 
(with transition rate 0). 



of our model, we could write equation ([7]) as, 

/■oo 

/ Aj(0,t\s ,0)dt= Aj(0,u = 0\s ,0) (10) 
Jo 



2/3 + 6 
2/3 + S + u 



2P + 5 + V 



Notice that as v grows (v >> 2/3 + 6) each Aj becomes 
smaller except for A\, with A\ — s> 1, i.e. the escape is 
certain in the first visit. However when v << 2/3 + 6, the 
probability of escape A(-) has contributions from each 
j-visit. This can best be understood considering, 



A j+ x(-) 



u=Q 



2/3 + 5 



U + 1) (11) 



which gives the relative contribution of successive terms 
in ([5]) . When v gets smaller the contribution is spread all 
over j values as Eq. (fir]) shows. In contrast, each Aj — > 
(J 7^ 1) as v grows, concentrating all the probability in 
A x . 



Mean Escape Time (MET) 

Following the ideas exposed in [l9| and by resorting to 
the matrix formalism and the Dyson's procedure f(26|) we 
were able to obtain the probability P(n, m, t\no, mo, t — 
0) (in the Fourier-Laplace space), which is the building 
block for the MET. For the detailed calculation see ap- 
pendix IA1 

We will denote the (finite) Fourier transform by its 
argument , as we did in the Laplace transform case. Thus 
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for example the transform on a coordinate, say x, would 
read: 

P(k, m, t\no, mo, 0) = J r {P(n, to, t\no, too, 0)} 

N-l 



n=0 



e ikn P(n,m,t\n o ,m o ,0) 



From P(k,m,u\nQ,mo,t = 0) (obtained in the Fourier- 
Laplace space), the probability that a walker is on the 
surface at site (n,m) at time t given it was at (no, too) 
at t = 0, P(n,m,t\no,mo,t — 0), is derived by us- 
ing the inverse Laplace transform on u and the inverse 
Fourier transform on k (for the x coordinate) for each 
[P(fc, u)] m . However, as we are interested in the calcu- 
lation of © , we only need to perform the inverse Fourier 
transform on P(0, 0, u\uq, mo, t = 0), i.e. we need the el- 
ements J- 1 |[P(fc,u)] mo |. We obtain for [P(fc, u)} mQ : 

\F(k u)] = 

' °' m ° 5(1 - T7)(l - Ty^- 1 ) + (M - Ax(fc))(l + 77^) ' 

(12) 

where ry = 1 + f u - \Ju 2 + 4-yuj /2j, M = 2M + 1 

and u = u — A{k). The inverse Fourier transform on 
[P(fc, u)] Q m is carried out in the following way, 



P(0,0,u|n o ,m o ,0) = — 

9 =o 



■[n 2 -^,u)] 0ma (is) 



Thus we have obtained the required expression for the 
calculation of the MET through the narrow escape win- 
dow and it only remains to choose the initial distribu- 
tion. We now evaluate the MET for a walker with an 
uniform initial distribution on the base line (y = 0). 
This means the initial distribution is given by g(n,m) = 
(1 — 5 n fi)S mi o/(N — 1). Notice that we explicitly exclude 
the possibility of having a walker at (0, 0) at t = [30| . 
This way we obtain, 



T=N 



M 

1 



N-l r 



TV- 



1 ^ 



2nq 
~N~' 



0) 



J 0.0 ) 

(14) 

We make some comments regarding equation (|14[) which 
constitutes one of our main results. Notice that (fT4|) 
adequately provides the limits of perfect escape case, v — > 
oo (obtained in [l9j]) and no escape window, v — > 0, T— > 
oo. Observe that as it is commented in [l6| for a perfect 
escape case, T could be expressed in the following way, 



T = 



NM 



7 
N 



.£„ [n 



2-rrq 
N ' 



°)]o,o 



N - 1 



2ft 



2ft 



E q [n 



2-rrq 
N ' 



°)]o,o 



N - 1 



(15) 



or T = T sur f ace +Ti ine ; (fT5|) provides an interesting physi- 
cal insight into the problem. Simply notice how the mean 



escape time is constructed from the mean duration of sur- 
face excursions and the mean duration of border or line 
excursions (first and second terms of (fT5|) respectively). 



Regarding the existence of a minimum in T 

T could be enhanced with respect to S provided we are 
able to find 5* — 6*(f3,~/,v, N, M) - the optimal desorp- 
tion probability - that satisfies, 



dT 
~85 



NM 
7^ 



N 



N 



^ Mj~ 1 2fta 2 -ai 
' N- 1 ( ai S* + a 2 2ft) 2 

q—1 x ' 



(16) 



where the a, = cti{q, N, M) for i = 1, 2 are defined in ap- 
pendix (jBj) . Notice that (fT())) defines an implicit equation 
for 5* which, although we could not solve, provided us 
with some generals conclusions as it approaches certain 
limits. Consider first ft — > 0, with the other parame- 
ters held fixed. In this case (a finite value for) S* exists 
whenever we have a finite escape probability rate v and, 



N-l 

\ M ^ N ~ J ) S ' 



(17) 



This is a highly interesting result since in the perfect 
escape case {y = oo) this minimum disappears, as 8* is 
pushed towards oo. Thus the 'imperfect' escape window 
enables a region that was absent in the perfect case. On 
the other hand in the limit ft >> 6 it could be shown 
that equation (|T6l) can not be satisfied for any 6* value. 
In this case, and taking into account the walker's initial 
distribution, the transport is performed on the baseline 
(lower boundary) of the confining domain, so this is an 
expected behavior. 

In the following section we will make more remarks 
regarding the minimum in T, while introducing some il- 
lustrations corresponding to S* . 



III. ILLUSTRATIONS 

Here we illustrate the general framework introduced in 
the previous section and compare our theoretical results 
to independent Monte Carlo simulations. In the next fig- 
ures, lines indicate analytical calculations while symbols 
correspond to Monte Carlo (MC) simulations. We will be 
interested in situations in which a mixed type of trans- 
port generates a global minimum in the Mean Escape 
Time. 

In Fig. |3]we present curves corresponding to the MET 
(Mean Escape Time), as a function of the desorption rate 
5, with parameters N — 20, M — 10 and ft = 0.1, for 
different values of the escape 'strength' v : which is the 
"transition rate" at the escape window. We have in- 
cluded for comparison the 'perfect escape case', i.e. once 
in the escape window the escape is instantaneous, with no 
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FIG. 3. MET as a function of the desorption rate (in log scale) 
S, with M = 10, N = 20, P = 0.1, for different values of the 
transition rate at the escape window v. From top to bottom 
v = 0.05,0.1,1,10,100. We have also included (thick solid 
line) for comparison, the perfect escape window case (y = oo). 



FIG. 4. MET as a function of the desorption rate (in log-log 
scale) 8, with /3 = 0.01, N = 10, M = 2, for different values 
of v (transition rate at the escape window) . From bottom to 
top v — oo, 100, 10, 1, 0.1, 0.05. Lines correspond to analytical 
calculations and symbols to Monte Carlo simulations. 



dwelling time. Notice how v regulates the existence of a 
minimum in the MET: as v gets smaller the imperfection 
rises the T curve until it becomes monotonous. Hence, 
for this situation we could say that v affects negatively 
the 'mixed' type of transport (journey's along the bound- 
aries and the surface). However it is worth remarking 
that the transition rate at the escape window can con- 
tribute positively as well. This behavior is well depicted 
in figure 2] 

Figure |4] presents curves corresponding to the Mean 
Escape Time, as a function of the desorption rate 8, with 
N = 10, M = 2, fi = 0.01, for different values of the tran- 
sition rate v. As can be inferred from the figure, v signifi- 
cantly influences MET as it varies from to oo. Changes 
in the location of the extrema values of MET can be seen 
ranging from a monotonous behavior (y — > oo extrema in 
S — >• oo) to a situation with a global minimum, and then 
back again into a monotonous behavior (y — > extrema 
in 8 — > 0). So in this case the transition rate to the es- 
cape window contributes 'positively' to the mixed type 
of transport, since it turns a situation without a mini- 
mum (perfect escape case) into a situation of enhanced 
transport (minimum in T for some values of v). 

In Fig(5]we present curves corresponding to the 8 value 
that minimizes MET , 8*, as a function of j3 for different 
values of i>, obtained from the numerical solution of Eq. 
(|16[) . All lines depicts quite a similar trend for finite v\ 
8*{p ,N,M, r ),v) = values marked by empty circles 
are not included in the curves and mark the end of the 



/3- interval in which 8* exists. In other words T is not 
monotonous while (3 S [0,/3 o ). As we show in Appendix 
(|B|) f3 satisfies, 



2/3 



7 



Z_=i aia 2 



M (JV-1)2&, 



(18) 



For values larger than j3 the T curve reaches a minimum 
at S = 0. However this is found at the beginning of 
the 5-interval and without change of sign of dT/88. We 
decided to rule it out as long as it doesn't represent a true 
interplay between surface and boundary paths. In these 
situations all particles stay in the base line and eventually 
escape trough the escape window without excursions into 
the surface. 

The behavior of 8* considerably changes in the perfect 
case (y = oo). Particularly the range of /3 values where a 
minimum exists in T shrinks as indicated by the dashed 
asymptotes in the figure. The left/right asymptotes indi- 
cate the limit in which T becomes monotonous, extrema 
for 8* — »■ oo and 8* — > respectively. The left and right 
asymptotes are respectively located at, 



N-i /N-i \ 

9=1 \q=l J 



(19) 
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FIG. 5. (color online) 8* as a function of /? (in lin-log;lin-log 
scale) for fixed N, M and 7 for different values of v. The 
cuts of 8* in the ft axis (8*(j3 ) = 0) for different «/, indi- 
cated by red empty circles are not included in the curve; fi 
points are obtained from equation (|18[) . Vertical dashed lines 
are the asymptotes for the perfect escape case obtained from 
equations |19| and (|20|) . 



For clarity's sake in the inset we have magnified the entry 
points to the j3 axis of S* curves for v — 10, 100, 00. 

Figure [6] shows the phase diagrams that summarize the 
existence /non-existence of enhanced transport, analyzed 
from the perspective of the existence of a minimum in 
the Mean Escape Time. The diagrams are plotted for 
fixed N, M and v as a function of the transition prob- 
ability over the baseline, /?, and the surface transition 
probability 7. White regions correspond to non-optimal 
transport (absence of minimum -monotonous behaviour- 
in the MET), while filled regions (red patterns) identify 
regimes with enhanced transport. We have also included 
in figure |51 enclosed by black lines the region correspond- 
ing to enhanced transport in the perfect window escape 
case. Dashed curves are obtained from the solution of 
equation [TBI with f3 > 0, 



2/3o 



-1 
2 



'E 



( ^EgQity 

2(N- 1) ^ \\' r M(N- 1) I 2(N-1) 



(21) 

As it was expected, when v grows the regions approach 
the perfect case and the escape, once in the window, be- 
comes certain and instantaneous. Notice that we obtain 
quite a good agreement between the region of optimal 
transport, evaluated from equation (|14[) . and the corre- 
sponding bounds derived from relations (I2ip for finite v 
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FIG. 6. (color online) Phase diagrams that summarize 
the existence /non-existence of enhanced transport for v = 
0.01, 1, 100, 00 respectively, for fixed system sizes N = 20 and 
M = 10. White regions correspond to non-optimal transport, 
while filled regions (red patterns) identify regimes of enhanced 
transport. Region enclosed by black lines correspond to en- 
hanced transport in the perfect trapping case (y — 00) while 
dashed (green) lines correspond to the bound -from eq, (|18|l - 
after which, T becomes monotonous. 



and OH), (EH) for v = 00. 



IV. CONCLUSIONS 

We have presented a model based on a master equation 
approach to the narrow escape time problem. In this 
study we introduced a finite transition probability, is, at 
the narrow escape window which allowed us to study the 
imperfect escape case. Varying from to 00, v allowed 
the study of both extremes of the trapping process: that 
of a highly deficient capture, and situations where escape 
is certain (perfect trapping case). 

By resorting to Dyson' technique we have obtained an- 
alytic results for the primary quantity studied in the NET 
problem, the Mean Escape Time (MET), and we have 
studied its dependence in terms of the transition (des- 
orption) probability over (from) the surface boundary, 
the confining domain dimensions, and the finite transi- 
tion probability at the escape window. Particularly we 
showed that the existence of a global minimum in the 
NET is controlled by the 'imperfection' of the escape pro- 
cess. Regarding such conclusion, a very interesting result 
was that the 'imperfect' escape window enabled a region 
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where T could be minimized, something the perfect case 
lacked. 



We have also presented bounds -equations (|T8|) . (|19|) 
and pop - between which an optimal minimum value of 
T could be found, improving previous bounds derived in 
[l9l | . The phase diagrams introduced in the last section 
deserve a special word, for not only do they present a 
compact summary of the situations of enhanced trans- 
port, whenever some exist, but they also can lead to a 
better understanding of the relations among the param- 
eters that characterize the system. In addition to our 
analytical approach, we have implemented Monte Carlo 
simulations, finding excellent agreement between the the- 
oretical results and simulations. 

We consider that the presented scheme is an analyt- 
ically manageable model, which could be used to study 
the impact of several (domain dimension, different rates 
of transition, etc.) parameters in the interplay between 
surface and boundary pathways, and could also serve as 
a forerunner for the study of more general and complex 
systems. This work contributes to an area of growing 
interest, providing a more general overview of a previous 
work (Pil) and showing a plausible physical insight into 
the surface-mediated diffusion mechanisms in the pres- 
ence of an imperfect escape window. 

The current approach to the narrow escape time prob- 
lem can be generalized in several directions: higher di- 
mensions, "dynamical" behaviour of the narrow escape 
window, non-markovian desorption, etc. All of these as- 
pects will be the subject of future work. 



(IH), we obtain: 

m = 

uP(k, 0, u) - P(k, 0,t = 0) = -yP(k, 1, u) 

-{8 - Ai(fc))P(fc,0,M) 
m = 1 

uP(k, 1, u) - P{k, l,t = 0) = 6P(k, 0, u) + jP(k, 2, u) 
-(2 7 - A(k))P(k,l,u) 
2 < to <M-1 
uP(k, to, u) - P(k, to, t = 0) = A(k)P(k, m, u) 

+7(P(fc, to + 1, u) + P(k, m-l,u) 
-2P(k,m,u)) 

m = M 

uP(k, M, u) - P(k, M, t = 0) = A(k)P(k, M, u) + 

jP(k, M - 1, u) - jP(k, M, u) . 

(Al) 

Here we have defined Ai(fc) = 2/3(cosfc — 1), A(k) = 
27(cos k — 1). Using the matrix formalism, equation (|A1[) 
can be written as 



[ul - H] P = I , 



(A2) 



where I is the identity matrix, H is an (M + 1) x (M+ 1) 
tri-diagonal matrix with elements: 



C; 7 

S C 7 ... 

7 C 7 

... '■• '■. '-. 







7 C 7 
7 7 + C* 



(A3) 
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C and C\ are defined as C — —2j + A(k), C\ = 
-5 + A 1 (k), and P is an (M + 1) x (M + 1) matrix with 
components, 



P(k, to, u\uq, mo, t = 0) 



P(M)] r 



In order to find the solution to equation (|A2j) we decom- 
pose the H matrix in the following way, 



= A(fc)i + H° + e 1 + e 2 , 



(A4) 



where 







"-7 


7 







" 


Appendix A: MET calculation 




7 


-27 


7 









e° = 





7 


-27 


7 





In this appendix we focus on the calculation of the 








7 


-27 


7 


probability P{n, m, t\n$, toq, t — 0), which is the building 











7 


-7. 



(A5) 



block for the Mean Escape Time. Taking the (finite) 
Fourier transform with respect to the x variable and the corresponds to the transition matrix for a symmetric ran- 
Laplace transform with respect to the time t in Equation dom walk to nearest neighbours in a finite lattice (M + 1 
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sites) with reflective boundary conditions at the ends. A formal solution to equation (IA2[) is: 

On the other hand: 

P = [ul - Hp 1 . (A8) 

i r a ^ applying the Dyson procedure ( [26jj ) a general ex- 
H = (T ~ 5 + ^i(fc)-A(fc))5 i>0 5oj=A 1 (5i,odo, i (A6) pre ssion -in the Fourier- Laplace space- for [P(fe, u)] m , mo 
H 2 = —(7 — <5) <5i 5 i 5o,j = A 2 5^1 Soj (A7) could be found, 
I 







^' ^ m,l + 1 _ A. M7i „M ( AQ ) 



l-(A 1 + A 2 )[P°(fc, u )] ,o \ L v ' l-AipP0(fc,«)] 0i0 



where. 



or 



[P°(M)] 



m,mo 



jy|m— too | _|_ ^M— (m+mo) 
27(1 — 77) -+- €t 



u = u - A(k), M = 2M + 1 and 77 = 1 + 

(u - \Ju 2 + A-yiij /2*y. 

From (jA9j) . the probability that a walker is at site 
(n, m) at time £ given it was at (noi TO o) at £ = 0, 
P(n, m, £|no, mo, £ = 0) is derived by using the inverse 
Laplace transform on u and the inverse Fourier trans- 
form on k (for the x coordinate) for each matrix element 
pP(*,«)] 



771, TOQ 



Notice that, as we are interested in the 
calculation of ([5]). we only need to perform the inverse 
Fourier transform on P(0, 0, u\tiq, mo, £ = 0) i.e., we need 

the elements J 7-1 | [P(fc, u)] rno j . In this case expression 

191) reduces to, 

yyTOo _^ „Af— too 

m ° ~<5(1 - tj)(1 - ry*- 1 ) + (u - Ai(fe))(l + ^) 

(A10) 

The inverse Fourier transform on [P(fe, u)] mQ is carried 
out in the following way, 



pP(M)] 0j 



JV-l 



P(0,0,u\n ,m ,t = 0) = jjYl 



9=0 



W— >«)]„ 



. (AH) 

Thus we have obtained all the required expressions for 
the calculation of the MET. We now proceed to eval- 
uate the Mean Escape Time for a walker with an uni- 
form initial distribution on the base line (y — 0), i.e. 
g(n, m) = (1 — 5 n ,o)5 m ,o/{N — 1). Notice that we explic- 
itly exclude the possibility of having a walker at (0, 0) at 
t = 0. We obtain, 



T = N 



M 

7 



v N - 



N-l 
9=1 



P(^ = 0) 



0,0 ) 
(A12) 



T = N 



M 1 

7 S 



N-l 



v N 



9=1 



1 



Q'i<5 + a 2 2/3 



where 



(A13) 



a 1 =a 1 ( q ,N,M)= i \^= o)i \j:= o) (A14) 



(l + ^=o) 2M+1 
2tt 



a 2 = a 2 (g, N) = 1 - cos —q 



(A15) 



Appendix B: 5* - optimal desorption probability 

In this section we present some results regarding the 
desorption probability rate value that minimizes T, 5*. 
For this recall equation (|16p . 



dT 
~85 



s=s* 



NM 



N 



N 



Mj- 1 2(3a 2 - ati 



= 



JV-l ^ (aid* + a 2 2/3) 2 

q=l x 7 

(Bl) 

Let us focus on the relation between 6 and /3 as these 
are the parameters of interest, since the former enables 
the transport on the surface, and the later regulates the 
movement on the boundary line where the escape window 
is located. Although we will keep track of all variables, it 
could be shown that (|B1|) can be written in terms of the 
scaled variables /?' = /?7 _1 , v' = v^ 1 and 5' — S-f^ 1 . 
So a modification in 7 would result in an enlargement or 
shrinkage (if we let 7 get smaller or larger respectively) 
of the former variables. We will consider the behaviour 
of d5* /dp. Even though we were not able to obtain an 
explicit expression for 8*, its derivate with respect to (3 
could be evaluated in closed form. To do this recall equa- 
tion (|B1I) and differentiate it with respect to j3 consider- 
ing S* = f(j3) (with the other parameters held fixed). 
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Then after some algebra we obtain, 



and rearranging some terms in (|B3|) we finally obtain 



OS* 
~dj3 



2(M~( 



S* + l)aia 2 



L^q=l 



(a 1( 5*+Q 2 2/3)2 



E 



N-l {M 1 - 1 &* + l)a 2 1 
3=1 (QK5*+a 2 2,3) 3 



■spiV-1 
Z^q=l 



(Q 1 A'*+a 2 2/3) 2 



(B2) 

From (|B2I) we could obtain the entry points to the j3 axis 
of 8* curves. The forerunner for this is the sharp growth 
on the 8* curves in figure ([5]). As a matter of fact the 
abrupt increase is in d\og8* /d\og (3. However is not dif- 
ficult to show that this happens only if the denominator 
in (|B2[) — > for some /3 (also notice that in this situation 
5* ~ 0) so, 



N-l 

E 



N-l 



=J (a 2 2/3 ) : 



= M7- 1 E 



a! 



=J (« 2 2/3o) 2 



(B3) 



-h 2 1 (iV — 1) 



JV-l 



2 ^°E 



where we have used (|B1[) to replace the sum on the right 
hand side. We could go even further and replace the 
sum on the left hand side. For this we differentiate 
(|B1|) with respect to 5* and get Mj~ 1 2(3J2 q atict2(ai5* + 
a 2 2/3)~ 3 = J2„ a i (a 1 8*+a 2 2/3)^ 3 . By using this relation 



2/?o = 



7 



E»=l "1^2 



M (JV-l)2/3o 



E»=l "2 



(B4) 



The solution of equation (|B4j) in terms of /3 that makes 
S*(j3 ,N,M,j,i') = marks the end of the interval in 
which 8* exists i.e., the T curve becomes monotonous. 

For the perfect escape case we obtain one of the asymp- 
totes between which 8* exits, outside them the extrema 
is pushed either to or to 00, by letting v — >• 00 in (|B4|) . 



2 As- 



7 2-,q=l a l a 2 

M Ef^" 1 



(B5) 



For the second asymptote we go back to (|B2|) . follow 
a similar reasoning that lead us to equation (|B3[) . here 
z/ = 00, 8* — > co, and obtain, 



T-^N-l -1 

7 Eg=i a i 



(B6) 



Equations (|B5[) and (IB6[) constitute an improvement to 
the bounds derived in [19| and the solution of equation 
(|B4[) gives a bound (we could not find similar results in 
the literature known by us) regarding the existence of a 
minimum in the MET. 
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